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Mathematical models of production of sheet glass were developed. The problem of controlling the melting- 
production problem was formulated. The efficacy of utilizing mathematical models for working out effective 
actions to improve the quality of the glass was demonstrated. 


The most important problem in the quality management 
system in production of sheet glass is constantly increasing 
the effectiveness of the manufacturing processes based on 
data analysis and elaborating the corresponding corrective 
actions. 

Simulated modeling of operation of the manufacturing 
line is being conducted to determine resources for further im¬ 
proving glass quality and different algorithms for controlling 
the glass production process are being investigated. Mo¬ 
deling was done with statistical data from operation of the 
1LPS line in 2005. 

Control of the sheet glass melting-production process is a 
multicriterial decision-making problem. The penalty func¬ 
tion method was used for solving it in [1], 

The selected solution is evaluated with a set of criteria 
which differ in relative importance: 

F = (Q % , ADe, B, Ri, Ze, Re, Tn, Cu, Bl, a), 

where Q is the gas consumption for glass melting, m 3 /h; 
ADe is the change in the density of the glass over 24 h, 
g/cm 3 ; B are the melting bubbles (class 1 and 2 defects), 
number per 10 m 2 ; Ri is rippling of the glass, evaluated by 
samples; Ze is waviness of the glass, measured with the 
Zebra method, deg; Re are the deviations in the reflected grid 
index, mm; Tn is the thickness nonuniformity of the glass, 
mm; Cu is the curvature (deviation from planeness), mm; 
Bl is the bloom effect, g/m 2 ; a are the residual stresses in the 
glass ribbon, nm/cm. 

Each local criterion characterizes some local goal of the 
decisions made. The optimum solution is selected in the re¬ 
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gion of compromises [2], Selection of the compromise 
scheme has not been formalized. The principle of selection 
of one optimized criterion was used in the study. The gas 
flow for glass melting, whose minimization in the given 
glass production plan corresponds to the problem of obtain¬ 
ing the highest gain in production, was selected with the opti¬ 
mized criterion. Restrictions are imposed on the other indi¬ 
cated criteria, and fulfilling them should ensure production 
of high-quality glass - M0 and Ml. 

Decision making based on correction of the glass pro¬ 
duction regime (control) is executed at input of the control 
object. The value of the density of the glass is established 
4 days after an action, while the optical distortions, parame¬ 
ters of the glass, and defects are determined during the cur¬ 
rent control step (24 h) due to the different inertness of the 
control channels based on regime variables. The lack of a 
strong correlation between the output variables (De, B, Ri, 
Ze, Re, Tn, Cu, Bl, a), allows combining planning of the 
density of the glass with prediction of the density several 
days in advance with the problem of controlling the optical 
properties and defects in the glass during the current deci¬ 
sion-making step. 

The glass quality indexes are described by linear regres¬ 
sion equations [1] in which the temperature of the gas me¬ 
dium in the furnace based on the readings of dome 
thermocouples © gm2 , © gm6 , © 8 along the axes 2nd, 6th, 
and 8th pair of burners; the temperature of the tin in the float 
bath bays © 12ba , © 20ba ; 1 ^ le temperature of the glass entering 
the float bath © en ; the temperature of the glass ribbon going 
out of the float bath © ex ; the speed of the first edge-forming 
machine v lefln ; the annealing temperature in annealing fur- 
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nace zones A and B - ® a and & h ; the slope of the tempera¬ 
ture curve along the annealing furnace tunnel a v 

The restrictions imposed on the properties and flaws of 
the sheet glass must be fulfilled: 

change in the density of the glass over 24 h: 

I De (i - 1) - De (z) | < ADe max , 

where De (0 gm6 , © gm8 , , C Si02 ); i = 2, 3, ... is the de¬ 

cision-making step with respect to control; 

the optical distortions visible in transient light: 

Ze ( v lefm > ®12ba ’ ©20ba ’ 5, C^) > Ze min ; 

the optical distortions visible in reflected light (grid): 

Re (®en > ©12ba > ©20ba’ ©ex ) - Re max 1 

the thickness nonuniformity of the glass ribbon: 

( v lefm ’ ®en ’ ©lba ' © 12ba ’ 5) — ^ n max ’ 

the deviation from planeness (curvature) of the glass 
sheet: 

Cu (§, 0 ex , De) < Cu max ; 
the Bloom effect: 

B1 (®i2ba > Cq 2 , © ex ) ^ Bl max ; 
the residual internal stresses in the glass: 

o(«,, v gr ,5)<a max , 

where Cj,^ and C' S|( ) are the content of Fe 2 0 3 and Si0 2 in 

the glass, %; § is the thickness of the glass ribbon, mm. 

It is necessary to select one of the many decisions on cor¬ 
recting the operating conditions of the equipment from the 
area defined by the system of restrictions: 

w gm2, min — w gm2 — w gm2, max’ 

©gm6, min “ ©gm6 “ ©gm6, max’ 

®gm8, min — ©gm8 — ©gm8, max’ 

©12ba, min — © 12ba — ®12ba, max’ 

©20ba, min — ®20ba - ©20ba, max’ 

0 < 0 < 0 

mm — a — a, max ’ 

©A, min — ©A — ©A, max ■ 

The magnitude of the increase in the controlling effects 
in each control step should not exceed the values determined 


by the features of the melting-production process techno¬ 
logy: 
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The restrictions on the shape of the temperature curves of 
the glass furnace and gloat bath should also be observed: 

(©gm6 0 ) ~~ ®gm5 ( z )) - A ©gm56’ 

(©gm6 0 ) ~~ ®gm8 (* )) - A ©gm68’ 

®12ba O' ) _ ©ex O') I = A ©12 e xi 

®20ba 0 ) - ©ex 0)1 = A ©20exi 

i ~ 1, 2, 3, 4,.... 

In solving the control (planning) problem, the sheet glass 
production plan and ribbon thickness are defined. The ribbon 
drawing rate v gr is determined by a calculation in terms of the 
planned production with consideration of the ribbon thick¬ 
ness and width. The regime variables not included in state¬ 
ment of the problem are stabilized at levels determined by 
the process regulations. 

Previously elaborated adaptive control models were used 
to solve the production control problem [1,3]. 

In consideration of the above, the following penalty 
function was obtained: 

M(t) = (X l abs (min (Re max - Re, 0)) + X 2 abs (min (Bl max - 
Bl, 0)) + A, 3 (abs (A0 2Oex - (©20ba (f ) — ©ex 0 )))) + 

(abs (A0 12ex - (0 12ba (t) - 0 ex (t)))) + 

(A 4 (abs (min(Aa - abs ( a(t ) - a(t- 1)), 0))) + 

*•5 ( abs ( min (( CT max + (5(0 - 2 )) - 0))) + 

A 6 abs (min (a, 20))) + (X 7 abs (min (ADe max - ADe, 0)) + 

hQgm + Z 9 abs ( min ( Ze - Ze min ’ °)) + 

*■10 abs (min (A0 gmfi) - abs (0 6gm (t - 1) - 0 6gm (O)> °)) + 

Oo abs (min (A0 gm8 - abs (0 8gm (f - 1) - © ggm (0)’ 0)) + 

Oo abs (min ((0 6gm (f) - 0 5gm (t) - A© gm65 ), 0)) + 

X w abs (max((0 8gm (O - 0 6gm (O + A0 gm86 ), 0))), 

where M(t) is the penalty function in the control step; A. is 
the penalty coefficient; t = 2, 3, ... is the decision-making 
step with respect to control. 
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TABLE 1 



Manual regime 

Control algorithm 


average value 

standard deviation 

average value 

standard deviation 

Density, g/cm 3 

2.4965 

0.0036 

2.5029 

0.0037 

Zebra, deg 

53.0 

6.0 

57.8 

6.0 

Gas consumption, m 3 /h 

5002.0 

362.7 

4821.0 

371.0 

Grid (1st section), mm 

4.28 

1.69 

3.25 

1.16 

Bloom effect, g/m 2 

1.250 

0.001 

0.890 

0.220 

Thickness nonuniformity (1st section), mm 

0.041 

0.033 

0.042 

0.025 

Residual stresses, nm/cm 

8.98 

2.00 

8.80 

3.40 

Curvature, mm 

0.026 

- 

0.009 

0.001 

Temperature of gas medium in glass furnace, °C: 
along axis of 6th pair of burners 

1511.5 

18.8 

1561.5 


along axis of 8th pair of burners 

1456.5 

17.7 

1410.4 

15.7 

Temperature of tin in float bath, °C: 
in 12 th bay 

798.8 

18.1 

813.2 


in 20th bay 

603.0 

5.0 

596.4 

- 

Temperature of glass ribbon coming out 
of gloat bath, °C 

616.0 

5.7 

612.5 

5.6 

Intensity of change in annealing temperature 
over tunnel length, °C/m 

-4.37 

1.7 

-4.34 

1.7 


Control of the melting-production process using the pen¬ 
alty function algorithm allows conducting the manufacturing 
process more economically in comparison to manual control: 
a decrease in consumption of natural gas by approximately 
4% is expected. The results of modeling the 1LPS line with 
2005 data are reported in Table 1. 

As simulation modeling showed, using the proposed 
control algorithm instead of manually conducting the melt¬ 
ing-production process makes it possible to more economi¬ 
cally utilize natural gas for melting glass while simulta¬ 
neously improving its quality. 


REFERENCES 

1. R. I. Makarov, E. R. Khorosheva, and S. A. Lukashin, Automa¬ 
tion of Sheet Glass Production (Float Method) [in Russian], Izd. 
Assotsiatsii Stroitel’nykh Vuzov, Moscow (2002). 

2. R. I. Makarov, V. V. Tarbeev, E. R. Khorosheva, et al., Sheet 
Glass Quality Control (Float Method): A Manual [in Russian], 
Izd. Assotsiatsii Stroitel’nykh Vuzov, Moscow (2004). 

3. A. V. Vasil’ev, “Development and study of mathematical models 
of sheet glass production units,” in: Proceedings of the Interna¬ 
tional Scientific and Practical Conference “Socioeconomic Sys¬ 
tems and Processes: Methods of Investigation and Problems in 
Development [in Russian], Vladimir (2005), pp. 325 - 328. 






